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Abstract. We extend our model of chemic al evo lution that successfully account for the main observables in the 
solar neighborhood (Alibes, Labay & Canal 2001) to the whole Milky Way halo and disk. We assume an inside-out 



scenario for the assembling of the Galaxy through the infall of external material. Two different compositions of 
the infalling material have been considered: primordial and slightly metal-enriched (0.1 Zq). Our calculations 
follow the evolution of 76 isotopes up to the Iron peak, but we only present here results for the elements having 
well-measured gradients: C, N, O, Ne, Mg, Al, Si, S, Ar and Fe. The calculated current radial distributions of 
these elements are compared with large samples of data from H II regions, B stars, planetary nebulae and open 
clusters. We also discuss the time evolution of the abundance gradients. Good fits are achieved with both infall 
compositions, closer to the observational data in the case of primordial infall. We show that contributions from 
intermediate-mass stars allow to reproduce the measured abundances of C and N, but the profiles of their ratios 
to oxygen leave still open the nucleosynthetic origin of these elements. In agreement with several previous works, 
our models predict a flattening with time of the abundance gradients, most of which taking place during the early 
galactic evolution. 
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'"i 1. Introduction 



Studies of the chemical evolution of the Galaxy have ben- 
efit in recent years from a growing number of high quality 
observations concerning chemical abundances in stars and 
in the interstellar medium (see Smartt 2001 for a review), 
from which important constraints on the ingredients of 
theoretical models can be inferred. In addition, on the 
theoretical side, new metallicity dependent yields derived 
from up to date stellar evolutionary calculations cover- 
ing the whole range of stellar masses have been provided 
by several groups (van den Hoek & Groc ncngc n 1997 ; 
Fores tini & Charbonnel \l99T\ Marigo et al. |1996| ; Marigo 
2001, for low and intermediate-mass stars; Woosley &; 



is assembled by continuous infall of primordial or lightly 
enriched gas, are able to reasonably reproduce most of the 
observational data in the solar ring (Boissier & Prantzos 
199$ Chang et al. |199$ Portinari & Chiosi |1999| , that in- 
cluded radial gas fl ows a nd the effe cts of a galactic bar in 
Porti nari & Chiosi 2000 ; Hou et al. |2000 ; Chiappini et al. 
|2001 ). Although different authors use different prescrip- 
tions for the basic ingredients (essentially the law for the 
star formation rate (SFR), the slope and the lower and 
upper mass limits of the initial mass function (IMF), the 
gas flows in and out of the Galaxy, and the nucleosynthetic 
yields) most of them share several basic features. In partic- 
ular, they usually consider an inside-out scenario for the 



Weaver 1995; Limongi et al. 2000, for massive stars). 

Among the different observables, the stellar metallicity 
distribution in the solar neighborhood (G-dwarf metallic- 
ity distribution) and the abundance gradients along the 
galactic disk play a major role. Standard closed-box mod- 
els of the chemical evolution of the galactic disk cannot 
adequately reproduce the G-dwarf metallicity distribu- 
tion, as they produced too many stars at low metallicities. 
However, the so called infall models, in which the Galaxy 
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formation of the galactic disk (Larson 1976) through radi- 
ally dependent timescales for both the infall and the SFR, 
IMF constant in space and time, and negligible galactic 
winds. 

However, the convergence reached by standard chem- 
ical models on the main observables in the solar neigh- 
borhood contrast with the results on the abundance gra- 
dients in the galactic disk. The existence of such gradi- 
ents, with single slope values dlog(X / H) / dr ~ —0.06 dex 
kpc -1 for most of the elements measured, seems today 
firmly established from observations of Hn regions, young 
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O and B stars, planetary nebulae (PN), and open clus- 
ters. Several simple models (not including radial flows) 
are able to reproduce the presently observed trend of the 
abundance distribution as a function of the galactocen- 
tric distance. Nevertheless, they substantially disagree re- 
garding the time evolution of the abundance gradients. 
Models by M olla e t al. ( |1997| ), Allen et al Q1998D , Boissier 
& Prantz os ( |1999| ), Portinari & Chiosi ( |1999[ ) or Hou et 
al. (200C) obtain steeper gradients in the past, while the 
opposite is found in papers by Tosi (1988) or Chiappini 
et al. (1997). This discrepancy results from the different 
radial dependence of the enrichment (SFR) and dilution 
(infall and mass loss from low-mass stars) processes of the 
interstellar medium adopted in the models. 

Unfortunately, the present observational data do not 
allow to discriminate between these opposite theoretical 
results on the history of the galactic abundance gradi- 
ents. In principle, measures of elemental abundances in PN 
could help to solve the discrepancy because different types 
of PN have stellar progenitors of different masses and ages. 
In the same direction, studies of abundances in open clus- 
ters, whose ages cover a quite wide interval, would provide 
clues to discriminate the history of the galactic gradients. 
There is some evidence, derived from studies of PN, for a 
steepening of the gradients in the last few Gyr (Maciel & 
Quireza 1999), but the uncertainties that still plague the 
ages attributed to those objects, the lack of good sam- 
ples in the outer Galaxy, and the difficulty to asses the 
importance of dynamical effects on the distribution of old 
galactic objects prevent us to extract any firm conclusion 
about the evolution of abundance gradients in the Galaxy 
from PN studies. On the other hand, iron abundances de- 
rived from open clusters of different ages point towards a 
fairly constant iron gradient for most of the Galaxy life- 
time, or even to a flattening of the iron profile during the 
last Gyrs. 

In this paper we address the chemical evolution in the 
Milky Way by extending our model for the solar neighbor- 
hood (Alibes, Labay & Canal |200lj hereafter ALC [200lD 
to the whole galactic halo and disk. We use a standard 
open model (radial flows are not considered), which treats 
the formation of the Galaxy by infall of external gas 
f ollow ing the double infall scenario of Chiappini et al. 
( |1997| ). Two different kinds of composition of the infall 
material are used: primordial and enriched. The enriched 
one is based on observations of High Velocity Clouds 
falling towards the galactic disk. Recent measurements of 
Wakker et al. ( |1999[ ) of a massive (~ 10 7 Mq) cloud falling 
into the disk of the Milky Way show that such objects 
are slightly metal-enriched (~ 0.09 Zq). In ALC |2001 



showed that we can successfully reproduce the main ob- 
servational constraints regarding the chemical evolution in 
the solar neighborhood: G-dwarf metallicity distribution, 
age-metallicity relation, current values of the star forma- 
tion and the supernovae rates, the gas and star surface 
density and the evolution of the main elements up to the 
Iron peak. As in our previous work on the solar ring, here 
we calculate the evolution of chemical abundances along 



the galactic disk for the 76 stable isotopes of the elements 
comprised between hydrogen and iron, but we will just 
show here the results obtained for those elements which 
have been extensively observed in the Milky Way. 

The plan of the paper is as follows. In Sect. 0we briefly 
remind the main ingredients of our model. In Sect. we 
present a short review of the currently available observa- 
tional data on the elemental abundance distributions in 
the Galaxy. Sect. compares with observations our re- 
sults for the two assumed compositions of the accreted 
matter. In Sect. 



4.1 



and Sect. 4.2 



we examine in detail the 
galactic radial profiles and compare our results with the 
most recent data on the galactic gradients derived fr om 
observations both of young and old objects. In Sect. 12 
the time evolution of the abundance gradients is discussed 
and compared with data from old objects. Finally in Sect. 
g| we summarize our main conclusions. 

2. Overview of the chemical evolution model 

As usual, we model the galactic disk by a succession of 
1 kpc wide concentric rings, extending from 2 to 20 kpc. 
Each ring evolves independently increasing its mass by 
infall of extragalactic material with primordial or slightly 
enriched chemical composition. No outflows nor radial 
flows are considered. Some combinations of infall, SFR 
and drift velocity pattern of radial flows can successfully 
reproduce the observed radial profiles in the galactic disk 
(Portinari & Chiosi 2000 ). Given our poor knowledge of 
the processes that induce those radial flows (viscosity of 
the gas layer, differences in angular momentum between 
the infalling gas and that in the gaseous disk, interactions 
with spiral density waves), their inclusion in our chemical 
models would increase the number of adjustable parame- 
ters. On the other hand, radial flows tend to steepen the 
chemical gradients, and the observed values set limita- 
tions to any radial mixing. Therefore, we prefer to remain 
within the frame of static models. We calculate the chemi- 
cal evolution of each ring by solving the standard integro- 
differential equations of chemical evolution (Tinsley 1980), 
relaxing the instantaneous recycling approximation. We 
use the same prescriptions for the infall law, the star for- 
mation rate (SFR), the initial mass function (IMF) and 
the chemical yields as in ALC2001. In the next subsec- 



tions we briefly describe them, mainly highlighting their 
radial-dependence characteristics . 

2.1. Infall 

Infall models of galactic chemical evolution ordinarily 
adopt a single exponentially decreasing law to describe the 
assembling of the galactic disk. Thus, this kind of models 
does not take into account the formation and evolution 
of the halo and thick disk components of the Galaxy, al- 
though the final results concerning the properties of the 
galactic disk only show small differences as compared with 
those obtained by more complex infall laws. In any case, 
for consistence with our previous work, we also adopt here 
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the two-infall model developed by Chiappini et al. (1997), 
that assumes that the Galaxy builds up through two suc- 
cessive accretion episodes. The first one is responsible for 
the rapid formation of the halo and thick disk, and is fol- 
lowed by a separate, second and more extended episode 
which assembles the thin disk on timescales that are ra- 
dially dependent, within the framework of an inside-out 
scenario for the formation of spiral disks. Hence, we adopt 
the following infall law: 



dcr(r,t) 
dt 



A(r)e 



-t/T T 



B{r)e 



-(t-4max)/TB(r) 



(1) 



where a(r,t) is the total surface mass density, tt is the 
timescale for the halo-thick disk phase, taken as 1 Gyr for 
all r, and td(t) is the corresponding timescale for the thin 
disk phase, that we assume to be a linear function of the 
galactocentric distance, according to the equation 



TB (r) = 0.923r - 0.846 Gyr 



(2) 



obtained by adopting an inner disk value of td (r = 2 kpc) 
= 1 Gyr, and requiring that in the solar ring (tq = 8.5 
kpc) td takes the value of 7 Gyr that gave the best results 
in ALC2001. Finally, t max , the time of maximum accretion 
into the thin disk, is set as 1 Gyr. 

The coefficients A(r) and B(r) in Eq. Ill are fixed by 
imposing that the current total surface mass densities of 
the thick and thin disks are well reproduced by the model: 

AM = ^tMg) 
V ' T T (1 - e-W-rr) 



B(r) 



[cr(r,t G ) -a T (r,t G )} 



r D (r)[l 



-(tG-*m a x)/TD(r)l 



where the current total surface density, tr(r, to), and that 
of the thick disk, cr^{r, t G ), are exponentially decreasing 
functions of the galactocentric distance with the same 
characteristic scale length of td = 3.0 kpc, and the present 
solar values are 54 and 10 Mq pc -2 , respectively. We set 
the Galaxy age as t G = 13 Gyr. 

As stated in the Introduction, we have calculated mod- 
els for two chemical compositions of the infalling mate- 
rial. In the first case, which we will term in the following 
"primordial" model, we consider accretion of primordial 
material during the whole lifetime of the Galaxy, while 
a second set of calculations (our "enriched" model) starts 
incorporating primordial matter during the quick phase of 
the halo-thick disk formation and then, in the thin disk 
epoch, they accrete material with a low (Z = 0.1 Zq) 
but nonzero metallicity, that being inspired by the obser- 
vations by Wakker et al. (1999) of high velocity clouds 
of low metallicity (Z ~ 0. 09 Zq ) falling onto the galac- 
tic disk. As shown in ALC2001, these two types of infall 



models produce comparable results for the characteristics 
of the solar neighborhood but, as we will see, they exhibit 
noticeable differences in the outer galaxy where, due to 
the long infall timescales, the star formation and the cor- 
responding metal enrichment of the interstellar medium 



have been quite low at all epochs, and metals in the in- 
corporated material tend to level off the radial abundance 
profiles. 

2.2. Star formation rate 

We adopt the same formulation for the SFR than in 
ALC[200l|. The global behaviour of the SFR depends on 



the surface gas density, <r g , according to the Schmidt 



( 1959| ) law, *& oc er™, with an efficiency coefficient that de- 
pends on the galactocentric distance through some power 
n of the local total mass surface density <r(r, t) 



*(r, t) = v 



' M pc" 2 Gyr" 1 



T n-\-m — 1 



(r@,t) 



(3) 



where cr(rQ,£) is introduced as a normalization factor in 
orde r to e xpress the efficiency constant v in Gyr -1 . As in 
ALC 2001 , we adopt a value of v = 1.2 which is the one 
that best reproduce the characteristics of the solar neigh- 
borhood. We take m = 5/3 and n = 1/3 since, according 
to Dopita and Ryder ( 1994 ), those values fit well the em- 
pirical link found by these authors in spiral disks between 
the H Q emission, tracing current star formation, and the 
/-band surface brightness, a measure of the contribution 
of the old stellar component and, therefore, of the total 
surface mas s den sity. We note that recent observations by 
Kcnnicutt (1998) of the correlation between average SFR 
and surface gas densities in spiral galaxies point towards 
an exponent in the Schmidt law of m ~ 1.5. It is also worth 
to mention that, as can be seen in Fig. nL when this SFR 
law is applied to the observed current gas (Dame 1993) 
and total mass (Rana 1991) density profiles, it precisely 
follows the observed pres ent S FR dis tribu tion across the 
galactic disk (Lyne et al. 1985 Rana 1991 ). 

An alternative treatment of the radial dep enden ce of 
the SFR ca n be found in Boissier & Prantzos ( 1999 ) and 
Hou et al. ( 2000|) , based on the theory of star formation 
induc ed by density waves in spiral galaxies (Wyse & Silk 
1989 ). Although we will only show here results obtained 
with the Dopita & Ryder law, we get very similar con- 
clusions from calculations that use the Wyse & Silk SFR 
formulation. 

2.3. Initial Mass Function 



Despite some claims of variations of the IMF in time 
and/or in place (see Eisenhauer |200l| ) , there are not yet 
solid enough reasons to abandon the common view of a 
universal IMF. Moreover, simple models of chemical evo- 
lution can only deal with average trends, and in this case 
constant IMFs are favoured, as noted by Chiappini et al. 
(2000J). Same as in ALC2001, in this work we pres ent rc - 



sults obtained using the IMF from Kroupa et al. ( |1993[) , 
which consists of a three-slop e pow er law. Compared with 
the classical IMF of Salpeter ( 1955|) , it is quite flat for the 
lowest masses (x — 0.3 for M < 0.5 Mq), it steepens to 
x = 1.2 for 0.5 Mq < M < 1 Mq, and for higher masses 
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Fig. 1. Current surface density profiles of total and gas 
mass, and the theoretical prescription for the SFR by 
Dopita & Ryder ( 1994 ) (normalized to the solar value) 
applied to these distributions, compared to observational 
estimates from Lyne et al. (1985) and Rana (1991) 



it has a x ~ 1.7 slope. We adopt a minimum stellar mass 
of 0.8 Mq and a maximum of 100 Mq, 

We have also calculated models that incorporate the 



IMF of Scalo (1986). The basic trends of our results are al- 
most unchanged, except for the fact that models that use 
the IMF for Scalo give slightly steeper abundance gradi- 
ents, but the differences are so minute that we will not 
discuss them further in this paper. 

2.4. Yields 

The contribution to the chemical enrichment of the inter- 
stellar medium by stars of different masses and metallic- 
ities, i.e. the stellar yields, is one of the most important 
ingredients of models of galactic chemical evolution. In the 
last years several groups have provided new yields for the 
complete range of stellar masses and for several metallici- 
ties ranging from to solar or even supersolar. The s tellar 



yields adopted in this paper are the same as in ALC pOOl 
We briefly outline them in the following. 

Yields for low- and intermediate-mass stars are from 



van den Hoek & Groenewegen (1997, hereafter HG1997). 
For massive stars we use the yields from Woosley & 
Weaver ( |1995 ; hereafter WW 1995)^1. Type la supernovae 
are includ ed according to the yields from Thielemann et 
al. ( |l993|) . We also include nova outburst nucleosynthe- 
sis (an important source of 13 C, 15 N and 17 0) through 
the work of Jose & Hernanz ( 1998J) . We note that, despite 
the efforts of the different groups working in this field to 
produce accurate stellar yields, all the currently available 
ones suffer still for severe uncertainties, linked to the poor 
knowledge of phenomena such as mass loss, rotation, con- 
vection, the reaction rate for 12 C(a,7) 16 0, modeling of 
core collapse supernovae, etc. 



We remind he re th at, as in our previous solar neighbor- 
hood model (AL Q2001 ), we take only half of the iron yields 
calculated by these authors. 



We also want to mention that, in the inner galaxy, high 
metallicities are reached and there are no published yields 
for stars of such high metal contents. The maximum values 



reached in the sources used here are Zn\ for the WW 1995 







yields, and 2 Zq in the case of those of HG1997. We con 



tinue to use these yields even when a galactic zone reaches 
higher metallicities, since we think that extrapolation of 
the stellar yields beyond these maximum values is a less 
reliable procedure. 

3. Observational data 

We concisely review in this section the main observables 
for the disk of the Milky Way, relevant to the chemical 
evolution, which we will take as for reference to compare 
our results with. As mentioned, for instance, by Boissier & 
Prantzos ( |1999| ), all the existing observational data only 
inform us on the present properties of the galactic disk, 
but not on its past history, excepting the temptative in- 
dications obtained from the study of abundances in PNII 
and PNIII (see below) or in open clusters. 

3.1. Gaseous, stellar and SFR radial profiles 

The current surface gas density distribution is taken from 
Dame ( 1993 ), who gives the observed surface profiles of 
atomic and molecular hydrogen in the Galaxy. The total 
gaseous profile is the sum of both, increased by 40% to take 
into account helium and heavier elements, and is shown 
by the dashed line in Fig. [j]. As is widely known, its most 
prominent feature is the molecular ring about 4-5 kpc 
from the galactic center. 

The observed stellar density profile is exponentially de- 
creasing outwards. Recent estimates of the characteristic 
scale length give values in the range 2.5 — 3 kpc (Sackett 

The stellar surface density at 
1 (Gilmore 



1997; Freudenreich 1998 



35 ± 5 M, 



pc 



the solar ring is a*(rQ,tG) 
et al. |198S| ). 

The SFR along the disk of the Milky Way, in terms 
of its local value, has been estimated from observations 



of the H a emissivity (Glisten & Mezger 1982 )^the distri- 
bution of supernova re mnant s (Guibert et al. 1978 ), that 
of molecular gas (Rana 1991 ) , or from the surface density 
of pulsars (Lyne et al. 1985 ), all of those methods giving 
results in good agreement with each other. We display in 
Fig. [y and g the data from the last two references. 

3.2. Abundance gradients 

Observational studies on the radial profiles of the abun- 
dances of several chemical elements have increased during 
the last decade. Most of them are related to young ob- 
jects like Hii regions and B-type stars, with ages clearly 
inferior to 1 Gyr. Thus, these objects only give informa- 
tion on the present day abundances in the galactic disk. 
Observations of PN and open clusters could, in principle, 
give us some insights on the past evolution of the abun- 
dance gradients. In Table [l] we directly collect most of the 
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data available in the literature on the radial gradients in 
the galactic disk derived from observations of H n regions, 
B stars and PN. Observed iron gradients are summarized 
in Table |. 

Since the pioneering work by Shaver et al. (f 983), the 
analysis of the emission lines from ionized regions around 
massive star associations, in the optical and the IR region 
of the spectrum, have been used to obtain chemical abun- 
dance determinations of O, N, S, and, in a lesser extent, 
of Ar and Ne, in H II regions located over a wide range of 
galactocentric distances (from to 18 kpc). The main un- 
certainties that affect these measurements are related to 
the determination of the electronic temperatures (derived 
by different techniques, varying from direct measurements 
using suitable diagnostic lines to photoionization models 
that require a priori assumptions for the flux of ionizing 
photons), and also to the importance of the temperature 
fluctuations in the gaseous nebulae, which have not been 
extensively considered. Most of these studies agree on a 
galactic gradient of O and S of ~ —0.07 dex kpc -1 , and a 
little bit steeper for N. Some authors (Vilchez & Esteban 



1996) suggest a flattening of the gradients in the outer 



Galaxy, but the most recent analyses favour a unique slope 
(Afnerbach et al. |l997j; Rudolph et al. J1997D- However, 



this general concordance has recently being challenged by 
Deharveng et al. (2000). They have determined O abun- 
dances in a sample of H II regions located at galactocentric 
distances between 6.6 and 18 kpc by means of a double 
temperature model of Hn regions. As in most of the pre- 
vious work, their data do not show any flattening of the 
O distribution in the external regions, but they obtain 
an oxygen gradient of —0.039 ± 0.005 dex kpc -1 , almost 
a factor of 2 natter than the tradition ally q uoted value 
but similar to that from Esteban et al. ( 1999] ). Therefore, 
even if the existence of abundance gradients derived from 
H II regions is unquestionable, more observational effort is 
needed to clearly establish their actual magnitude. 

During the early 90's, the measurements from Hn re- 
gions showed much steeper gradients than those obtained 
from the photospheric composition of young B-type stars, 
which are supposed to reflect the composition of the mate- 
rial from which they formed, except for the rapidly rotat- 
ing B stars, where there could have been some mixing of 
atmospheric and core material (Maeder 1987). Spectra of 
those bright stars show optical lines of C, N, O, Mg, Si and, 
in high-quality spectra, also Al, S, and Fe. A series of stud- 
ies (Gehren et al. 1985; Fitzsimmons et al. 1992); Kaufer 



et al. |1994| Kilian-Montenbruck et al. |1994| ) detected no 



gradients or gradients much flatter than those determined 
from H II regions. Nevertheless, recent analyses of homo- 
geneous samples of B stars in a consistent manner (a cru- 
cial factor to obtain reliable measures when using stars) 
by Smartt & Rolleston ( |1997[ ), Gummers bach e t al ( |1998D , 
Rolleston et al. ( 2000[ ), and Smartt et al. ( 2001 ) confirmed 
the existence of gradients similar to those measured in 
i onize d nebulae. Furthermore, the study by Smartt et al. 
( 2001 ), which included four stars near the Galactic Center, 
found little evidences of O abundance variations inside 9 



kpc, i.e. a flattening of the O profile in the inner disk, at 
odds with the rest of the studied elements. 

The analysis of emission lines in planetary nebulae 
(PN) is another popular method to measure abundances in 
different regions of the galactic disk. Besides their greater 
statistical importance owing to the larger samples avail- 
able as compared with studies of H n regions and B-type 
stars, PN could give, in principle, information about the 
evolution of the abundance profiles. Following the orig- 



inal classification of PN by Peimbert (197S), revised by 
Pasquali & Perinotto (1993), disk PN are classified into 
three types according to their chemical and dynamical 
characteristics. PNI are He and N rich, belong to the thin 
disk, show peculiar velocities below 60 km s , and are 
associated with 2.4 - 8 Mq progenitors, thus having ages 
less than 1 Gyr. PNII are also members of the thin disk 
with low peculiar velocities, but they are not particularly 
enriched in He nor N; it is supposed that they come from 
stars in the mass range 1.2 - 2.4 Mq, which correspond 
to ages comprised between 1 and 7 Gyr. As for PNIII, lo- 
cated in the thick disk with velocities above 60 km s -1 , 
their assumed progenitors are stars in the mass range 1 - 
1.2 Mq, these then being the oldest disk PN, with ages in 
excess of 7 Gyr. However, large uncertainties are involved 
in the precise determination of distances, ages, and masses 
of PN. In addition, there exists the possibility that the ra- 
dial distribution of old disk PN might have been altered 
by dynamical effects through diffusion processes along the 
disk that could affect the gradients determined from those 
objects. Also, it is not clear whether PN show exactly the 
metal abundances of the gas at the time of their birth. 
Some enrichment by the products from the CN cycle affect 
all types of PN. In the case of PNI, whose progenitors are 
massive intermediate-mass stars, a possible enrichment 
by ON cycle products may exist, which could induce mis- 
leading oxygen abundance measurements in th ese nebu- 
lae. Rec ent st udies (Pasquali & Perino tto 1993 ; Maciel & 
Koppen 1994 ; Maciel & Quireza 199S ) determined gradi- 
ents of O and S about —0.06 dex kpc -1 for the available 
sample of PNII (representative of the thin disk) and a bit 
flatter, ~ —0.04 dex kpc -1 , for Ne. Some indications of 
shallower gradients in the outer (r > 10 kpc) Galaxy exist 
in those data, as well as a slight tendency of the gradients 
to flatten with time, but due to the considerable uncer- 
tainties noted above none of these points can be taken as 
definitive. 

Finally, abundances of Fe and its galactic gradient have 
been mainly determined by means of spectroscopic and 
photometric indices of stellar populations in open clus- 
ters. Early determinations (Janes 1979); Cameron 1985) 



found gradients of —0.06 to —0.09 dex kpc over a range 
of galactocentric distances from 7 to 16 kpc, a trend con- 



firmed by Friel & Janes (1993). However, Twarog et al. 
( 1997 ) , in an analysis of a sample of 63 clusters spanning 
a wide interval of ages, did not find any linear trend for 
the Fe abundances in the galactic disk. These authors sug- 
gest that a step function, instead of a unique slope, would 
better describe the distribution, with the discontinuity lo- 



(i 
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Table 1. Calculated abundance gradients in the galactic disk, in dex kpc 1 , and observed values obtained from Hn 
regions, B stars and planetary nebulae 



Ar 



C 



N 



O 



Ne 



Mg 



Al 



Si 



Hn 



Shaver et a. 1983| 



Ar 



Simpson et al. 



1995 



Vilchez & Esteban 



1996 



Afflerbach et al. [1997 
Rudolph et al 
Esteban et al. 



1997 



199E 



Deharveng et al. 2000 



5-14 
0-10 
12-18 
0-12 
0-17 
6-9 
6-18 



-0.133 

±0.022 



-0.090 
±0.015 
-0.100 
±0.020 
-0.009 
±0.020 
-0.072 
±0.006 
-0.111 
±0.012 
-0.048 
±0.017 



-0.070 
±0.015 



-0.036 
±0.020 
-0.064 
±0.009 



-0.049 
±0.017 
-0.039 
±0.005 



-0.080 
±0.020 



-0.045 
±0.017 



-0.010 


-0.060 


±0.020 


±0.015 


-0.070 




±0.020 




-0.041 




±0.020 




-0.063 




±0.006 




-0.079 




±0.009 




-0.055 


-0.044 


±0.017 


±0.030 



B stars 



Gehren et al. 1985 



Fitzsimmons et al. 
Kaufer et al. 11994 



1992 



1994 



Kilian-M. et al 
Smartt & Rolleston 
Gummersbach et al 
Rolleston et al 



1997 



1998 



2000 



Smartt et al. 2001 



8-18 
6-13 
6-17 
6-15 
6-18 
5-14 
6-18 
0-18 



-0.026 

±0.009 

-0.001 -0.017 

±0.015 ±0.020 



-0.035 
±0.014 
-0.070 
±0.020 
-0.070 
±0.020 



-0.078 
±0.023 
-0.090 
±0.010 
-0.060 
±0.020 



-0.010 
±0.020 
-0.030 
±0.020 
-0.000 
±0.009 
-0.021 
±0.012 
-0.070 
±0.010 
-0.067 
±0.024 
-0.067 
±0.008 



-0.043 -0.020 -0.002 ±0.000 -0.026 
±0.011 ±0.011 ±0.019 ±0.018 ±0.025 



-0.082 -0.045 -0.107 

±0.026 ±0.023 ±0.028 

-0.070 -0.050 -0.060 

±0.010 ±0.015 ±0.010 

-0.090 -0.050 -0.060 

±0.020 ±0.010 ±0.010 



PN 



1994 



1994 



Pasquali & Perinotto 1993 

(PN I - II) 
Maciel & Koppen 

(PN I) 
Maciel & Koppen 

(PN II) 
Maciel & Koppen 

(PN III) 
Maciel & Quireza 

(PN II) 



1994 



199£ 



Martins & Veigas |200C 
(PN II) 



1-16 
4-13 
4-13 
4-13 
3-14 
5-13 



-0.050 
±0.010 



-0.084 
±0.034 



-0.030 
±0.010 
-0.030 
±0.007 
-0.069 
±0.006 
-0.058 
±0.008 
-0.058 
±0.007 
-0.054 
±0.019 



-0.050 
±0.020 
-0.004 
±0.008 
-0.056 
±0.007 
-0.041 
±0.008 
-0.036 
±0.010 
-0.069 
±0.034 



-0.075 


-0.060 


±0.008 


±0.008 


-0.067 


-0.051 


±0.006 


±0.006 


-0.063 


-0.034 


±0.010 


±0.010 


-0.077 


-0.051 


±0.011 


±0.010 


-0.064 




±0.035 





Model results 


Enriched model 


4-16 


-0.051 


-0.066 


-0.047 


-0.043 


-0.052 


-0.061 


-0.052 


-0.054 


-0.057 




4-10 


-0.045 


-0.054 


-0.040 


-0.035 


-0.040 


-0.041 


-0.039 


-0.040 


-0.039 




10-16 


-0.064 


-0.085 


-0.062 


-0.058 


-0.070 


-0.089 


-0.072 


-0.076 


-0.082 


Primordial model 


4-16 


-0.058 


-0.074 


-0.053 


-0.048 


-0.061 


-0.077 


-0.058 


-0.062 


-0.065 




4-10 


-0.049 


-0.058 


-0.043 


-0.038 


-0.045 


-0.047 


-0.043 


-0.045 


-0.044 




10-16 


-0.076 


-0.100 


-0.071 


-0.065 


-0.085 


-0.117 


-0.082 


-0.088 


-0.094 



cated at 10 kpc. Further studies by Carraro et al. (1998) 
did not shown evidence of any abrupt discontinuity and 
their data fit a unique slope relation, in agreement with 
earlier measurements. A similar conclusion was reached 



by Friel (1999), using high-resolution abundance deter- 
minations for metallicity calibrators, who finds a slope of 
~ —0.06 ± 0.01 dex kpc -1 , shallower but still compatible 
with previous values (Table S summarizes the Fe gradi- 
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Table 2. Calculated iron gradients and observations in 
open clusters 



Open Clusters 



Ac 



Fe 



Janes 197£ 



198 



Cameron 

Friel & Janes 11993 



Piatti et al. 1995 



Twarog et al. 1997 



Carraro et al. 1998 



Friel 199S 



8-14 


-0.050 ±0.008 


7-11 


-0.110 ±0.020 


8-16 


-0.095 ±0.017 


7-14 


-0.070 ±0.010 


6-16 


-0.067 ±0.008 


6-10 


-0.023 ±0.017 


10-16 


-0.004 ±0.018 


7-16 


-0.085 ±0.008 


7-16 


-0.060 ±0.010 



Model results 



Enriched model 


4-16 


-0.059 




4-10 


-0.055 




10-16 


-0.069 


Primordial model 


4-16 


-0.064 




4-10 


-0.059 




10-16 


-0.079 



ents derived from available data). In any case, our current 
knowledge on the iron gradient as derived from open clus- 
ters is far from being clear. When comparing these results 
and those corresponding to extreme Pop I objects (Hn 
and B-type stars), we must note that open clusters can 
be as old as 8 Gyr, so that they do not trace the young 
component of the galactic disk. It is well known that Fe is 
mainly a product of Type la supernovae and thus its abun- 
dance history is very different from that of other elements, 
like oxygen and sulphur, which are almost exclusively syn- 
thesized in massive stars. The fact that the slope of the Fe 
gradient is similar to that of O is therefore surprising, and 
might be simply coincidental. To make things even more 
confused, the slope of the Fe gra dient in open clusters 
shows no dependence on age (Friel 1999 ), in other words, 
its gradient has not evolved significantly during most of 
the lifetime of the disk. This fact is in clear contradiction 
with the well-known age-metallicity relation from F stars 
(Edvardsson et al. |l993t Rocha-Pinto et al. ^000[ ). When 
in the next Section we compare the calculated Fe abun- 
dances across the disk with measures in open clusters, 
the reader should remember all the uncertainties affecting 
these data. 



4. Model results versus observational constraints 



We showed in ALC2001 that our models nicely reproduce 



the chemical evolution of the solar neighborhood for both 
compositions of the accreted matter, although a little bet- 
ter fit was obtained for the enriched composition. In this 
section we compare our results for the primordial and 
the enriched infall models with the available observational 
data on the properties of the Milky Way disk. 



4.1. Radial profiles: gas, stars and SFR 

The results obtained for the radial distributions of gas, 
stars and SFR are displayed in Fig. g. Our model reason- 
ably follows the Dame (1993) data on gas surface den- 
sity, <jg, although, as it usually happens in models that do 
not consider radial flows nor the dynamical influence of a 



galactic bar (Portinari & Chiosi 1999, 2000), the peak in 
gas density corresponding to the molecular ring at 4 kpc is 
clearly less marked than the observed one. The calculated 
stellar profile, which is exponentially decreasing, is also in 
very good correspondence with the observations. We ob- 
tain a radial scale length of 2.8 kpc, just in the range of 
2.5-3 kpc given by Sackett ( |1997[ ). Finally, the observed 
current star formation rate normalized to its solar value 
(\I//\I/q) is nicely reproduced by our calculations outside 
4-5 kpc, while in the most internal part of the galactic 
disk, even if still compatible with the observational limits, 
we obtain values a bit lower than the data. 

In summary, the model gives a good representation of 
the main characteristics of the galactic disk, in particu- 
lar outside 4 kpc, which corresponds to the regions where 
most of the data on the abundance gradients have been 
obtained. Therefore, we can confidently go on with the 
discussion on the metal abundances across the disk. 

4.2. Current abundance gradients and profiles of 
abundance ratios 

We begin by discussing the present abundance gradients 
and the profiles of abundance ratios to oxygen for the 
chemical elements with measured abundances. Their time 
evolution will be presented in Sect. 4.3. In Fig. || to || we 
show the radial distributions of C, N, O, Ne, Mg, Al, Si, S, 
Ar and Fe, as well as the profiles of their abundance ratio 
to O, obtained for both the enriched (solid lines) and the 
primordial (dashed lines) models, and we compare them 
with most of the existing data from young objects (H II 
regions, early B-type stars and PNI). In the case of Fe, its 
abundance has essentially been measured in open clusters, 
thus data in the Fe graph cover a wide interval of ages. 
The calculated values of the abundance gradients for all 
the chemical elements considered appear at the bottom 
of Table [[], except for Fe which are included in Table |[ 
We give the gradients over an interval of galactocentric 
distances of 4 to 16 kpc, thus covering both the inner 
Galaxy (4-10 kpc) and the outer Galaxy (10-16 kpc). 

Models of galactic chemical evolution that include a 
radial dependence of the SFR and the infall rate (the 
"inside-out" scenario), like ours, naturally produce more 
or less steep abundance gradients across the disk. We re- 
mind that the parameters of our model were tuned to 
fit the solar neighborhood constraints, in particular the 
chemical abundances of the Sun at the moment of its birth 



4.5 Gyr ago (see ALC j2001|) . It is well known that the solar 
neighborhood is more oxygen-rich than the local interstel- 
lar medium. This is the reason why the majority of the 
calculated radial abundance distributions shown in Fig. |3J 
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Fig. 2. Results for the gas and star surface density profiles, 
and normalized current star formation rate distribution 
compared with observations. Data: gas surface density 



from Dame (1993) accounting for 40% of He and the ob- 
servational error; star surface de nsity using <7*(rQ,tG) = 
35 ± 5 M Q p er 1 ( Gilmore et al. |l989|) and r* = 2.5 - 3 
kpc (Sackett |l997|) ; current SFR from Lyne e t al. (|l985|) , 
from the distribution of pulsars, and Rana ( 1991 ), from 
the distribution of molecular gas 



the Galaxy, the metallicity is modest even at the present 
epoch, and models with enriched infall give larger abun- 
dances of metals and shallow gradients as compared with 
those obtained when the accreted matter is assumed to be 
of primordial composition. 

The inner disk abundance profiles are flatter than in 
the external regions, a common result in most models 



of galactic chemical evolution (Portinari & Chiosi 1999 



Boissier & Prantzos 1999; Hou et al. 2000). The reason for 
this flattening has already been discussed in those papers: 
due to the initially high SFR in the central regions, the 
metal enrichment is appreciably diluted by the low metal- 
licity gas ejected by the large amount of long-lived, low- 
mass stars formed at early times from metal-poor gas, an 
effect that can not be compensated by the ejecta of mas- 
sive stars at late times, when very little gas is left and the 
SFR is drastically reduced. This is not the situation in the 
outer disk, where the surface gas density remains low dur- 
ing the whole evolution and the population of old stars is 
scarce. It must be noted that the central flattening of the 
abundance profiles is overestimated by our models. The 
WW 199?; stellar yields for massive stars only cover up to 
solar metallicity, and the HG1997 yields for intermediate- 
mass stars just arrive to twice solar. Then, the use of these 
yields when metallicities higher than solar are reached in 
the central zones underestimates the metal enrichment, 
especially for those elements that have stellar yields that 
depend on the stellar metal contents (odd-even elements, 
for instance). 

Since most of the studied elements are produced in 
massive stars and have stellar yields hardly dependent on 
Z, i.e. they are primary elements, we get flat [X/O] ratios 
across the galactic disk. The two exceptions are Al, whose 
yields clearly show the odd-even effect, and N which is 
mildly secondary in our intermediate-mass stars. 

In the remaining of this subsection we will briefly 
discuss our results for each individual clement. 



to gjonly fit the more metal-rich measurements. Although 
several suggestions have been put forward to account for 
the solar overabundances (orbital diffusion, dust deple- 
tion, enrichment of the protosolar nebula by neighboring 
SN II...), there is no convincing in contradiction with pre- 
dictions from standard chemical models where metallicity 
is an increasing function of time. However, whatever the 
real oxygen abundance in the present solar ring may be, 
what specially counts for our discussion is the radial be- 
haviour of the abundances, i.e. their gradients, since the 
particular election of the zero-point does not alter the ra- 
dial characteristics. 

We obtain a general good agreement with the data. 
Results for both compositions of the infalling matter are 
quite similar in the inner disk, because the high metallic- 
ity reached there, due to the high efficiency of star for- 
mation, erases any possible influence of a non-primordial 
composition of the accreted matter. However, in the outer 
disk, where the SFR remains low during the whole life of 



Oxygen.- Because O is a good metallicity indicator its 
radial abundance profile in the Milky Way disk, as well 
as in external disk galaxies, has been extensively stud- 
ied. Despite of early discrepancies between results from 
galactic Hn regions and PN, which give steep gradients, 
and observations of B stars, that showed much flatter 
gradients, there seems nowadays to be a general conver- 
gence towards an abundance gradient with a unique slope 

dlog(0/H)/dr 0.07 ± 0.01 dex kpc" 1 (see Table |). 

However, there are some hints that the situation might be 
less clear. The O gradient in the inner Galaxy as traced 
by young B stars appears to be rather flat (Smartt |2001 ) 
in contradiction with results from H n regions (Afflerbach 
et al. 1997). The reasons for this discrepancy are unclear, 



specially so in view that the other elements considered in 
the same paper do not show such inner flattening. A fur- 
ther conflict comes from the recent work of Deharveng et 
al. ( 2000| ) concerning Hn regions (an indicator that tra- 
ditionally has given steep gradients). These authors have 
collected a new sample of Hn regions located at galac- 
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Fig. 3. Calculated current O abundance profile in the 
Milky Way disk. Solid line: enriched infall model. Dashed 
line: primordi al inf all model. Dat a: cir cles-H.il regions 
(Shav er et al. 1983 ; Fich fc Silk ey 1991 ; Simpson et al 



1995; Vilchez & Estcban 1996; Afflerbach et al. 1997 



Rudolph 1997 ; Dehar veng et al. 200C ); squares-H stars 



1998 



(Smartt & R ollcsto n 1997 ; Gumme rsbac h et al 
Hibbins et al. |l998| Rolleston et al. ^000| Smartt et al. 
2001 ). Sol ar Sy stem abundances (©) are from Anders & 
Grevesse Ql989|) 



tocentric distances in the range 6.6 to 17.7 kpc, and de- 
termined O/H abundances assuming a two-temperature 
model for the ionized nebula. They obtain a much shal- 
lower O gradient (by almost a factor of two) than previous 
surveys of Hn regions. 

Our models reasonably reproduce the O abundance 
profile as displayed in Fig. |3J. There is an small zero-point 
discrepancy, as stated above, and we get slightly shallower 
gradients than those derived from B stars. The enriched 
model gives a result in good agreement with the value 
quoted in Deharveng et al. ( [2000 ). As it is the case for 
the other elements, the primordial infall model yields 
a steeper slope, which is still a bit smoother than the 
canonical value, but compatible with observations, in 
particular those of young PN. In both models the O 
profile in the inner disk is flatter than in the external 
zones. 

Carbon.- Measures of C abundances in young disk ob- 
jects are few and rather uncertain. Observations of Hn 
nebulae in a small range of galactocentric distances (6-9 
kpc) by Esteban et al. (199E) give a steep carbon radial 
gradient of -0.133 ± 0.022 dex kpc -1 . Most of the addi- 



tional data come from B stars. Gummersbach et al. (1998) 
found a much smoother gradient of —0.045 ± 0.014 dex 
kpc -1 . However, l ater s tudies by Rolleston et al. (2000) 



and Smartt et al. ( 2001 ) extending from the inner galaxy 
to 18 kpc gave gradients of the order of —0.06 ± 0.02 
dex kpc -1 . Results for the C abundance distribution and 
the corresponding [C/O] profile are presented in Fig. [|. 
Our enriched model produces d\og(C/H)/dr ~ —0.05 dex 
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Fig. 4. Present C abundance and C abundance ratio to O 
(referred to the solar value) profiles. Solid line: enriched 
infall model. Dashed line: primordial infall model. Data: 



squares-B stars (Gummersbach et al. 1998; Hibbins et al. 
1998; Smartt et al. 2001 ). Solar System abundances (©) 
are from Anders & Grevesse ( 198S| ) 



kpc -1 , while in the calculation with a primordial compo- 
sition of the accreted matter we obtain a steeper value of 
~ —0.06 dex kpc -1 , in excellent accord with the recent 
data sets. The local carbon abundance at £q = 8.5 Gyr 
is nearly solar, which shows that contributions from low- 
and intermediate-mass stars are enough to complement 
the carbon produced by the WW 1995 yields for massive 
stars. 

Nevertheless, our models produce a flat C/O profile 
(see the lower panel of Fig. |I|) because carbon production 
hardly depends on metallicity (or even decreases as the 
metallicity grows) in the HG |1997| yields. Measures of this 
ratio in local low-mass stars point to an increase of C/O 
with metallicity (Gustafsson et al. |1999 ). In fact, in the 
recent study of a sample of B stars by Smartt et al. ( 200l[ ) 
a gradient of —0.05 ± 0.02 for C/O is found, steeper than 
the —0.02 ± 0.03 previously measured by the same group 
(Rolleston et al. 2000] ), owing to the inclusion of four 
galactic center stars, carbon-rich but not oxygen-rich, 
in the sample of Rolleston et al. (2000). In this last 
paper, the authors themselves notice that their absolute 
[C/O] ratios are significantly offset from the Orion and 
solar values. They attribute this discrepancy to errors 
in LTE C and O abundance determinations that do not 
invalidate the underlying trend of their [C/O] ratios. If 
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Fig. 6. Present Mg, Si and S abundance and abundance ratio to O profil es. S olid line: enriched infall model. Dashe d 
line: primordial infall model. Data: circles-Hn regions (Affl erbac h et al. 1997 ); squares-B stars (Smartt et al. 2001 ) 
Solar System abundances (0) are from Anders & Grevesse 



this trend is real, since it cannot be explained by the 
contribution to carbon abundances by intermediate-mass 
stars, we must turn to massive stellar yields that include 



infall model), in good agreement with the observational 
determinations. Low- and intermediate-mass stars con- 
tributions are mandatory to match the local abundance, 



mass loss, such of those of Maeder (1992), where C is 
produced at a higher rate than O in metal-rich stars 
due to the enhanced mass-loss in the post He-burning 
phase. Recent papers by Henry et al. (2000]) and by 



as pointed out by Hou et al. (2000) who, including only 
yields from massive stars, found subsolar N abundance 
in the solar ring, even when the evolution was calculated 



Hou et al. (200C) have shown that metallicity-enhanccd 
carbon yields from massive stars d la Maeder ( 1992 ) 
could account for the observed carbon abundances in the 
galactic disk, producing at the same time C/O profiles 
increasing towards the galactic center. In any case, to 
settle the question of the role played by stars of different 
mass ranges in carbon enrichment, more extensive and 
accurate determinations of the C/O ratio across the 
galactic disk are needed. 



Nitrogen.- N has been extensively observed in Hll re- 
gions, B stars and PNII (see Table |]). In general, its gra- 
dient seems steeper (~ —0.08 dex kpc -1 ) than that of O, 
but there is a wide scatter in the values quoted in the lit- 
erature, going from —0.05 dex kpc -1 , as determined by 
Pasquali & Perinotto (|1993|) in PN, to -0.11 dex kpc" 1 



with the Maeder (1992) stellar yields 



in the study of Hn regions by Rudolph et al. (1997). The 
recent analysis by Smartt et al. (2001), including N abun- 
dance determinations in four young B-type stars lying 
toward s the Galactic Center in the sample of Rolleston 
et al. ( |2000|) , derived a N gradient of -0.06 ± 0.02 dex 
kpc -1 . In our calculations we also obtain a steeper N gra- 
dient as compared with the O profile, with values ~ —0.07 
dex kpc -1 (as always, slightly steeper in the primordial 



The ratio N/O is essential to understand the relative 
importance of primary versus secondary production 
of N. A clear tendency of N/O ratios increasing with 
metallicity (i.e. with decreasing galactocentric radius) 
would indicate that secondary production dominates 
de N evolution. Observations of external spiral galaxies 
(Vila-Costas & Edmunds |l993j ; Henry et al. |2000| ) show 
that at metallicities higher than 12 + log(0/H) ~ 8.3 
nitrogen behaves as secondary. In the Galaxy, H II regions 
data clearly show higher than solar N/O ratios in the 
inner Galaxy, but nearly solar values in the outer disk 
(Vilchez & Esteban |1996| ; Rudolph et al. |1997[ ). However, 
observations of B stars showed a unique and rather steep 
slope, from -0.04 ±0.02 dex kpc -1 ( Rolles ton et al. |2000D 
to -0.06 dex kpc -1 (Smartt et al. |2001[ ), for the N/O 
radial profile. Contribution to nitrogen abundances by 
low- and intermediate-mass stars allows, in our models, 



to rightly fit the solar value. The N yields of H GJ1997 
present a moderate secondary behaviour, that translates 
into a more pronounced N gradient than in the case of 
O. Therefore, our models give N/O profiles increasing 
towards the inner Galaxy, characterized by a gradient in 
N/O of ~ —0.02 dex kpc -1 , as shown in the lower panel 
of Fig. Q Nevertheless, this result is still much less than 
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Fig. 7. Present Al, Mg and Ar abundance and abundance ratio to O pro files. Solid line: enrich ed inf all model. Dashed 
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the values derived from B stars. Thus, our calculations 
support the idea that intermediate-mass stars give an 
important contribution to the building up of the nitrogen 
abundance in the Galaxy, but new measurements of 
the N/O ratio are crucial to asses whether the nitrogen 
synthesis in those stars has a clear secondary character. 



The a elements Magnesium, Silicon and Sulphur. - 

These three a-elements should correlate with O since all 
of them are produced and returned to the ISM through 
Type II supernova explosions. The observed radial profiles 
of Mg and Si have been obtained from B stars surveys, and 
they are well fitted by linear relationships with slopes of 

0.07±0.01 and 0.06±0.01 dex kpc" 1 , respectively, 

compatible with the O gradient at 1 a level. S has been 
observed in Hn regions and PN, and its gradient is also 
similar to that of O. Our results show a nice agreement 
with the data (see upper panels of Fig. || and Table [l|), in 
particular in the case of our primordial infall model, which 
again gives steeper gradients than the enriched model. 

We display our results for the [X/O] ratios for the 
three elements in the lower panels of Fig. ^. As expected 
on theoretical grounds, their calculated abundance ratios 
to O along the galactic disk are flat, with solar absolute 
values for Si and S, while in the case of the Mg/O ratio we 
obtain slightly subsolar values in the whole disk, alt hough 



the difference is not large. As pointed out in ALC2001 



and also by Goswami & Prantzos ( |2000D , the WWJ1995 
yields underproduce Mg, showing, besides, unexplained 



metallicity effects. 



Aluminum, Neon and Argon.- Our results for these 
elements are presented in Fig. R (upper panels: gradi- 
ents; lower panels: [X/O] ratios). Aluminum has been ex- 
clusively measured in B- type stars. Gummersbach et al . 
( |1998|) , Rolleston et al. (|2000|) and Smartt et al. ( |200l| ) 
found relatively shallow Al gradients, of the order of 
—0.05 ± 0.01 dex kpc -1 . This is rather surprising, since 
Al is an odd monoatomic element ( 27 A1), subjected to the 
well known odd-even effect, i.e. stellar yields increasing 
with metallicity. Therefore, owing to the low metal enrich- 
ment of the outer galaxy, the Al abundances there should 
increase more slowly than those of pure primary elements 
like O, giving rise to a rather steep gradient. In fact, that 
is precisely the behaviour that we obtain in our models 
(see Fig. ^). Owing to the marked odd-even effect in the 
Al yields of WW 1995, we produce an Al gradient steeper 
than those of O and the a-elements. As a consequence, 
although our calculated present local Al/O ratio is almost 
solar, it shows a noticeable dependence on galactocentric 
distance, and hence metallicity, in the outer Galaxy (more 
evident in the case of our primordial infall model) , again 
at odds with the data that show no trend at all. Similar 
results, with even a more pron ounced Al gradient, were 
obta ined by Hou et al. ( |2000 ), who also employed the 
WW199J yields for massive stars. In view of the above, 
one is forced to conclude that theory and observations do 
not match each other, and that the current Al yields for 
massive stars should be revised. 
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Fig. 5. Present N abundance and N/O profiles. Solid 
line: enriched infall model. Dashed line: primordial infall 



model. Data: circles-Mil regions (Shaver et al. 1983; Fich 
fe Si lkey |199lt Simpso n et a l. |1995| ; V flchcz & Esteban 
1996| ; Afflerbach et al. |1997| Ru dolph |1997[ ); squa res-B 
stars (Gum mersb ach et al. 199? ; Rolleston et al. 2000; 



Smartt et al. 2001 ). Solar S ystem abundances (©) are from 
Anders & Grevesse (|l989 ) 



Ne and Ar have been mainly measured in PN, and 
in much lesser extent in Hn regions. The large set of 
data from PN (Maciel & Koppen 1994; Maciel & Quireza 
1999) give a Ne gradient of ~ —0.04 dex kpc -1 , and 
a somewhat larger value of ~ —0.05 dex kpc -1 for the 
Ar gradient. Our models produce results for these two 
elements in very good agreement with the observations, 
and comparable to those calculated for the rest of the 
a-elements. The corresponding ratios to O are ~ solar 
in the local ring, with the expected flat distribution. 
Nevertheless, the Ne/O profile remains below most of 
the observations. Regarding this last point, we remind 
that PNI are probably O-poor, and that measurements 
of Nc /O in PNI could be overestimated (Peimbert et al. 
19951 ). 



Iron.- In Fig. || we present the calculated Fe profile com- 
pared to the current data, while the values of the cor- 
responding slopes appear in Table |2|. Fe abundances are 
usually measured in open clusters. Earlier determinations 
of its radial gradient gave values in the range —0.06 to 
—0.09 dex kpc -1 . Although there have been claims that 
a two-step distribution with an abrupt discontinuity at a 
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Fig. 8. Calculated iron profile compared with observations 
in open clusters. Solid line: enriched infall mo del. D ashed 
l ine: p rimordial infall mod el. Da ta: Cameron (1985 ), Friel 
( |1995| ) and Carraro et el. ( |l998| ) 



galactocentric distance of ~ 10 k pc would best adjust the 
obser vation s (Twarog et al. 1997 ), recent work by Carraro 
et al. ( 199S|) and Friel (199£) support linear radial Fe pro- 
files, with a slope similar to that found in the former stud- 
ies. We obtain a Fe gradient of ~ —0.06 dex kpc -1 , in ex- 



cellent agreement with the value quoted by Friel ( 1999 ) . It 
is disputable to what extent open clusters trace the current 
abundances profiles, since those objects span a rather wide 
interval of ages. Nevertheless, when studying open clusters 
samples of different ages there is a general agreement on 
no or just very little evolution of the Fe gradient during 
the galactic disk lifetime. Actually, however, the absence 
of any age-metallicity relation for the galactic open clus- 
ters, in clear contradiction with the local field stars, casts 
some doubts on our understanding of the Fe evolution in 
the Galaxy. 

4.3. Gradient evolution 

A crucial point relative to the gradients of elemental abun- 
dances in the Galaxy, that can help to discriminate be- 
tween theoretical models, is their time evolution. Have 
they been steepening or flattening? Or we just have had 
a constant gradient during most of the Galaxy evolution? 
As indicated in the Introduction, most models that suc- 
cessfully account for the majority of the chemical features 
of the solar neighborhood, as well as of the present galactic 
disk, find a gradual flattening with time of the radial abun- 
dance distributions (e.g. Portinari & Chiosi 1999; Hou et 
al. 2000). Nonetheless, there also exist models that give 



abundance gradients that progre ssivel y steepen with time. 
For instance, Chiappini et al. ( 2001 ), with a model in- 
volving star formation thresholds for both halo and disk, 
obtain abundance gradients that increase in the course of 
galactic evolution. Among other reasons, this opposition 
between predictions from different models comes basically 
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from the diversity in the prescriptions used for the time 
dependence of the SFR and the infall. 

Unfortunately, the present status of observations can 
not help very much in solving the controversy. We have al- 
ready mentioned that abundance measurements obtained 
from PN or open clusters could give some insights on 
this subject, owing to the relatively wide range of ages 
covered by those objects. Maciel & Quireza (1999) found 
that abundance gradients derived from PN are shallower 
by 0.01-0.03 dex kpc -1 than those obtained from young 
objects, and roughly estimated a slowly steepening rate 
of —0.004 dex kpc -1 Gyr -1 for the oxygen gradient, but 
large uncertainties are involved in this result. First, the 
ages, masses and distances of PN are crudely determined. 
Second, abundances in PN might not be truly representa- 
tive of the composition of the ISM from which their parent 
stars formed. Finally, the distribution of PN across the 
disk could have been altered by orbital diffusion. 

The same way as the age-metallicity relation for 
nearby stars is a record of the metal enrichment in the 
local Galaxy, Fe measures in open clusters with a variety 
of ages do reveal some of the past history of the galactic 
abundance profiles. Besides, it is expected that orbital dif- 
fusion hardly affects the open cluster population, and their 
ages are more reliably determined than those of PN. Most 
studies indicate that the mean iron cluster abundance 
shows no dependence on age, pointing to a prompt initial 
enrichment of the disk (Friel 1999). This trend, i.e. that 
the radial iron gradient has not or just slightly changed 
with time, is found as well in Carraro et al. (1998). In fact, 
this last work hints to somewhat steeper gradients in the 
past, for intermediate-age clusters, although due to the 
small sample analyzed this result has little statistical sig- 
nificance. Results in the same directi on ha ve been prelim- 
inarily obtained by Bragaglia et al. ( 1999J ), though again 
with poor statistics. In conclusion, open clusters temp- 
tativcly indicate that the radial gradients have remained 
constant, or even slightly flattened, during most of the 
disk lifetime. 

We display in Fig. || the radial distribution of the ten 
elements that we are discussing, obtained with the primor- 
dial infall model, at different times. We do not present our 
results for the enriched model since they are quite similar, 
apart from giving, as mentioned above, shallower gradients 
in the outer Galaxy. The gradient evolution is more clearly 
seen in Fig. QG, where the explicit time evolution of the 
abundance gradients in the 4-16 kpc region is shown, both 
for the enriched (solid line) and the primordial (dashed 
line) models. We see a rapid flattening of the gradients for 
the majority of the elements, with the main changes in the 
radial gradients being produced during the first two Gyrs. 
Obviously, this trend is most noticeably when we consider 
enriched infall, since at early times the metallicity of the 
matter that settles into the galactic disk quickly erases 
the initial gradients. The evolution is almost completed 
at 5 Gyr. Afterwards, the flattening proceeds at a very 
slow pace till reaching the current values listed in Tables 
[l] and g. We must keep in mind that the enriched infall 
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Fig. 9. Abundance profiles for different times (1, 5, 8.5 
and 13 Gyr) in our primordial infall model. Long dashed 
lines correspond to the time of formation of the Sun. Solar 
system abundances (0) from Anders & Grevesse ( 198E ) 



is a kind of upper limit for the metallicity of the accreted 
material. A value of 0.1 Zq is what is observed nowadays, 
but the infall could have been metal-poorer in the past. 
In such a case, the gradient evolution would be something 
in-between those displayed for the two extreme compo- 
sitions considered. The abundance profiles in the inner 
Galaxy remain flatter during most of the evolution, due 
to the initially high infall and star formation rates, which 
induce a prompt saturation of the ISM, whose metallicity 
practically does not increase in later epochs (besides the 
fact that the stellar yields that we use do stop at solar 
metallicity, as already mentioned). 

The aluminum gradient evolution in the primordial in- 
fall model is clearly different from the rest, in the sense 
that the gradient steadily decreases during the whole evo- 
lution. The reason is once again the aluminum nature by 
being an odd-even element, with yields that increase with 
metallicity. This effect is seen as a slower evolution of the 
gradient as a consequence of the idle enrichment in the 
outer Galaxy. 

In summary, our models, as it is the case in some previ- 
ous work on this topic, produce abundance gradients that 
tend to decrease in the course of the Galaxy lifetime, ba- 
sically in its earlier stages. This result is in concordance 
with the indications about the Fe abundance behaviour as 
deduced from studies of open clusters, but it is in contra- 
diction with evidences obtained from PN. Actually, the 
current observational status on the temporal evolution 
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Fig. 10. Evolution of the gradients in the range 4-16 kpc. 
Solid line: enriched infall model. Dashed line: primordial 
infall model. A quick flattening of the gradient is obtained 
in the enriched infall model before 2 Gyr, due to the accre- 
tion of 0.1 Zq material. At £q it still produces smoother 
gradients than the primordial infall model 



of the chemical abundances, both in the Galaxy and in 
other spirals, is subjected to such uncertainties that it is 
far from allowing us to discriminate whether the galac- 
tic abundance profiles steepen or flatten with time. For 
instance, our model predictions can reasonably accommo- 
date the data from PN. Following the work of Hou et al. 
(200C), in Fig. Owe compare, in the case of the primordial 
infall model, the abundance profiles at various epochs with 
the observed distribution in PNII (ages between 1 and 7 
Gyr) and PNIII (ages > 7 Gyr), for the four elements (O, 
Ne, S and Ar) that have been measured there. We reach 
a mild agreement with the data from PNII for Ne and Ar 
and in less extent for O, while S appears clearly above the 
data. In any case, the spread of PNII abundances, that 
could in part be attributed to the errors associated with 
distance and age determinations, is a little bit larger than 
the range of our predictions. The PNIII abundances, how- 
ever, are well reproduced by our model, as it is apparent 
in the right panels of Fig. [0], where almost all the data 
points fall between our predicted abundance profiles at 
times of 6 Gyr (age of 7 Gyr) and 0.1 Gyr (age of 12.9 
Gyr). 



5. Summary and conclusions 




Fig. 11. Abundance profiles obtained in our primordial 
infall model at different times compared with data from 
PN. Left panels correspond to PNII data (with ages be- 
tween 1 and 7 Gyr) and the abundance profiles are from 
our models at t = 6 and 12 Gyr. Right panels are for 
PNIII observations (ages > 7 Gyr) and model results for 
t = 0.1 and 6 Gyr. PN abundances from Maciel & Koppen 



The chemical model developed in ALC [200l| , that nicely 
accounts for the observational characteristics in the solar 



(1994) and Solar System abundances (©) from Anders & 
Grevessc ( |l989| ) 



neighborhood, has been extended to the whole galactic 
disk. Our model assumes an inside-out scenario for the 
assembling of the Milky Way disk through infall of exter- 
nal gas, with a timescale that is a linear function of the 
galactocentric distance. Two different compositions of the 
accreted material have been considered. In our primor- 
dial infall model we assume that the matter that settles 
down to build the disk has a composition corresponding to 
the Big Bang nucleosynthesis. We have also investigated 
models, our enriched infall models, in which the infalling 
material that forms the thin disk has been preenriched up 
to a mctallicity of 0.1 Zq, in line with metallicity mea- 
surements of High Velocity Clouds. The adopted star for- 
mation rate has an explicit dependence on the total sur- 
face mass density, so that it implicitly becomes a function 
of the radial distance. Using metallicity-dependent yields 
for the whole interval of stellar masses, we have calculated 
the evolution of all the stable isotopes between hydrogen 
and zinc, although in the present paper we only analyze 
the elements that have been observed in the Galaxy, for 
which we have collected from the literature a wide and re- 
cent set of data. We have concentrated on the study of the 
radially-dependent properties of the galactic disk, giving 
special emphasis to the abundance profiles. Our models do 
not include radial flows because we prefer to keep as short 
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as possible the number of adjustable parameters (but see 
Chiosi 



Portinari 



2000 



for a recent discussion on the ef- 
fects of radial flows in standard chemical models). 

We briefly resume our principal conclusions as follows: 



We obtain satisfactory fits to the observed gas, star 
and SFR distributions in the Milky Way disk, although 
we do not reproduce adequately the molecular ring at 
4 kpc, a common feature of models that do not take 
into account radial flows and/or the effects of a central 
bar. Nevertheless, we fairly adjust the abundances of 
gas and the level of star formation in the galactic disk 
beyond the molecular ring, precisely the region where 
most of the elemental abundances have been measured 
and, consequently, the crucial zone for determining the 
chemical gradients. 

All the studied elements show noticeable gradients at 
the present epoch, whose values are compatible with 
the data from young objects, like H n regions and B- 
type stars, at lu level. For the best observed element, 
oxygen, our calculated gradient is a little bit flatter 
than the measured one, but still within the observa- 
tional limits. Models that include enriched infall give 
shallower gradients than models that assume a primor- 



— Our model, as it also happens in other work that 
adopts similar assumptions, predicts that abundance 
gradients in the disk do flatten in the course of galac- 
tic evolution. If the material that settles into the disk 
has an enriched composition, most of this flattening oc- 
curs during the firsts 2 Gyr, while this period is more 
extended when we consider primordial infall. Later on, 
the gradients continue to flatten but in a quite gentle 
manner. In other words, they exhibit little evolution in 
recent times. Secondary elements depart to some ex- 
tent from this behaviour, in the sense that their gradi- 
ents do not hastily smooth out, and they show a more 
steady pace towards their present values. 
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